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Abstract 

A two-dimensional integrable system being a deformation of the rational Calogero- 
Moser system is constructed via the symplectic reduction, performed with respect to 
the Sklyanin algebra action. We explicitly resolve the respective classical equations 
of motion via the projection method and quantize the system. 
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1 Introduction 



The relationship between the integrable systems and the Lie group gauge symmetries is 
well known However, Lie groups by no means exhaust all the symmetries and 

in fact there are a lot of systems whose gauge symmetries cannot be reduced to the Lie 
group action . For this reason it is quite interesting to construct an integrable system 
using the non Lie group symmetries. The aim of the paper is to present the example 
of such integrable system. We obtain the system via the symplectic reduction from an 
extended phase space, where the symmetries are realized by the action of the quadratic 
Sklyanin algebra Q]. Let us briefly explain the key idea for realizing this action. Note 
that the well known coadjoint action of a Lie algebra on its dual space can be written in 
terms of the linear Poisson structure, naturally associated to the Lie algebra. Namely if 
Pfj, are coordinates in the dual space and the respective Poisson tensor is of the form 
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where f^^ are the structure constants of the Lie algebra, then the coadjoint action is 
rewritten as 

SePt^ = atJiu{pY{p), (1) 

where are the infinitesimal parameters of the action. If we assume now that the Poisson 
tensor a^i,{p) in is no longer linear in we get the generalization of the coadjoint 
action to the case of an arbitrary Poisson manifold. Such a generalization of the coadjoint 
action has been originally proposed by Karasev ||^ and in general it cannot be reduced to 
the Lie algebra action in a sense that one cannot choose the global basis of the sections 
e^{p) with a Lie algebra commutation relations of the transformations (|1|). Instead the 
form of the commutation relations in a general case looks as follows 



[6^, 6rj] = 6, 

where 

e = V] = da{^, T]) + a{d^, rj) + dr]), 

(2) 

^ = 

dPf, 

Thus we get an example of the manifold, equipped with the non Lie algebra symmetries. 
A mathematical rigor treatment of such symmetries involves a definition of the so-called 
Lie algebroid 0,0, which is defined as a bundle A over a manifold M with a Lie bracket, 
given in the set of sections r(^) and an anchor map 6 : r(^) ^ r(TM) , being a 
homomorphism of the respective Lie algebras. Using this terminology we see that the 
above construction is an example of the Lie algebroid. Namely, the Lie algebroid is the 
cotangent bundle over the Poisson manifold. The Lie bracket between the sections C, and 
7] is given by the equation (|^) and the anchor map is defined by the formula (|l|)0. Note 
that the Jacobi identity of the Lie bracket holds in virtue of the respective Jacobi 
identity for the Poisson tensor a^,y(p) . In this paper we start from the Poisson manifold 
R"^, equipped with the standard classical Sklyanin bracket , and thus equipped with 
the algebroid action, which may be regarded as the action of the Sklyanin algebra. We 
choose the extended phase space of the integrable system to be the cotangent bundle 
over with the standard symplectic structure. We lift the algebroid transformations 
from the base manifold to the canonical ones in the cotangent bundle and readily find 
the respective Hamiltonian generators, which define a natural analogue of the momentum 
map for the case of the non Lie group symmetries. Then we construct the phase space 
of the integrable system by performing the symplectic reduction on the surface of the 
non zero level of the generators. The underlying Hamiltonians are obtained by reduction 
of the Casimir functions of the Sklyanin bracket. Note that since the Sklyanin bracket 
is the deformation of the Poisson bracket, which is equivalent to the linear bracket, 
associated to the Lie algebra U2 the presented symplectic reduction is similar to the 
one leading to the rational Calogero-Moser system |[l[],[jrT|. In particular, one of the 



"such example of the Lie algebroid, naturally associated to a Poisson manifold has been originally 
presented in the paper ||l0|| . 
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exponential degenerations of the Sklyanin bracket leads to the system, whose dynamics is 
in a sense equivalent to the dynamics of the two-particle rational Calogero-Moser system. 
As an aside, it is worth noting that although the integrable system is obtained via the 
symplectic reduction we still cannot find the respective Lax representation because the 
symmetries, that have been used for constructing the system are not Lie group ones. 
However, we show that it is possible to resolve the equations of motion for our system 
via the projection method. So the problem of finding the respective Lax representation 
is seemed to be rather intriguing and we are going to comment on it in the concluding 
section. The organization of the paper is as follows. In the second section we present 
the construction of the integrable system. In the third section we develop the projection 
method and explicitly resolve the equations of motion. The fourth section is devoted 
to the quantization. We find that although it is not possible to pose the bound states 
problem for our system one of the Hamiltonians possesses a discrete spectrum. At last, 
in the concluding section we discuss some open questions. All through the paper the 
Greek indices take values from to 3 the Latin indices run from 1 to 3 and whenever the 

indices i,j,k are meat they are assumed to take values 1,2,3 together with their cyclic 

d d 
permutations, d'^ denotes — — and is reserved for 7- — . 

2 Construction of the System 

We construct the system starting with the extended phase space, which is chosen to be 
the cotangent bundle T*R^ , equipped with the canonical symplectic structure 

uj = dpf,A dx^. (3) 

Here we use the following unusual conventions. The coordinates on the base are 
denoted by and the coordinates in the cotangent space are denoted by x'^ . Let us 
endow the base manifold with the quadratic Sklyanin bracket 

[Po^Pi] = o:oiip) = JjkPjPk {Pi,Pj} = Oiijip) = PoPk, (4) 
where the coefficients J12, J23 and J31 satisfy the identity 

J12 + J23 + ^31 = 0, (5) 

which is equivalent to the Jacobi identity for the bracket (^). The condition (|^) also 
implies that these coefficients can be represented in the form 

Jij Ji Jjy (6) 

where the values Jj are defined up to the additional constant c, Ji 1-^ Jj + c . Following 
the steps, outlined in the introduction we firstly define the analogue of the coadjoint 
transformations on the base manifold R^ using the Poisson structure (^. Namely, the 
infinitesimal coadjoint transformations are defined as 



SeP^ = aM'^(p)e''(p), 



(7) 
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where e'^ are the respective infinitesimal parameters. The transformations (|^ can be 
easily lifted to the canonical ones in the cotangent bundle 

SeP^. = «M-e"(p) S,x^ = d^{a,xe''{p)x^) (8) 

and the Hamiltonian generators corresponding to the transformations (|]) can be also 
readily found[| 

M, = a{e,x) = a^^ip)e''xr (9) 

Then in order to obtain a non-trivial dynamical system via the symplectic reduction 
method we consider the surface of non-zero level of the generators (Q) in the extended 
phase space T*R^. In analogy to what it is done for case of the rational Calogero-Moser 
system we set 

Mo = J23P2P3X'^ + J^lP^PlX^ + Jl2PlP2X^ = 
Ml = Po{x'^P3 - X^P2) - J23P2P3X^ = ^ 

(10) 

M2 = Po{x^Pi - xVs) - JsiPsPix^ = 

M3 = Po{x^P2 - x'^pi) - Jl2PlP2X^ = 0. 

The symplectic form being reduced on the constraint surface ( ]I0[ ) is of course degen- 
erate and the respective kernel distribution is given by the following vector field, being 
tangent to the surface (|TDp 

Spo = J23P2P3 Sx^ = (x'^ps - X^P2) 

5pi = 5x^=0 

(11) 

Sp2 = -P0P3 Sx'^ = -{poX^ + J23P3X^)e^ 
Sp3 = P0P2 Sx^ = {pox'^ - J23P2X^). 

Our system is then presented as a result of the symplectic reduction on the constraint 
surface (0) with respect to the transformations, corresponding to the vector field (|ll]). 
In what follows, these transformations are called the gauge ones and the integral curves 
of the vector field ( |TT|) are called the gauge orbits. We perform the symplectic reduction 
using the gauge fixing condition x^ = 0. In fact it is easy to see that this gauge condition 
has a horizon, given by the equation 

Pox^ - J32P3X'' = . (12) 



Namely, for the points, satisfying (12) the coordinate is unchange under the gauge 
transformations (|l3); and hence for these points the value of the coordinate x^ in general 
cannot be reduced to zero with the help of the gauge transformations. However, in the 

''thus we have constructed the so-cahed Hamiltonian algebroid over r*R'' ; in the paper jl^ it is 
shown that the natural construction of the Hamiltonian algebroid exists for an arbitrary Lie algebroid 
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paper we do not bother the subtlety and we are not going to discuss it. So taking into 
account the condition = together with the constraint equations (|1^) we get the 
embedding of the reduced phase space into the extended one. Namely, the embedding is 
given by the following mapping 

= = = = 

n ^ (13) 

where u^,Vi,u^, are the canonical coordinates parameterizing the reduced phase space. 
The pair of the Casimir functions of the Sklyanin bracket 

C, = \ Y.(Pif C2 = \{{po? + E Mpi?) (14) 
I ^ I 

provides us with the required set of the commuting and independent Hamiltonians of 
the presented system. Really, the gauge invariance of the Casimir functions and their 
involution property with respect to the canonical structure on T*R^ ensures that the 
Hamiltonians, obtained by reduction of (|1^ on the surface of the constraint equations 



and of the gauge fixing condition are commuting and independent. Thus we readily find 
the Hamiltonians for the our system 

- ^ + 2[V,U- - h2u'v2Y + ^^^^ 

where we set J2 = using the arbitrariness of the coefficients Jj . We have presented 
the example of the integrable system, obtained via the symplectic reduction, which is 
performed in the context of the non Lie group symmetries. The most similar symplectic 
reduction performed in the group theoretical context is the one leading to the two-particle 



rational Calogero-Moser system |jl|,|]TT|. This resemblance of the reductions becomes quite 
clear if we compare the Poisson manifold (^) to the dual space of the Lie algebra U2 
and the coadjoint transformations (|^ to the coadjoint action of the Lie algebra U2 . Note 
that the exponential degeneration of the Sklyanin bracket corresponding to case J23 = 
leads us to the system with the following Hamiltonians 

which describe the dynamics of the two-particle Calogero-Moser system with the coupling 
constant, depending on the initial value of V\ and with v} , being the coordinate of the 
mass center. 



3 Projection Method 

In the this section we show that in spite of the more complicated nature of the gauge 
transformations the equations of motion for the presented system ( plSf ) can be explicitly 
resolved via the projection method. Let us briefly recall the main idea of the method. It is 
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well known that each point of the reduced phase space corresponds to a gauge orbit on the 
constraint surface, where the dynamics is described by the gauge invariant Hamiltonians 
(([l^ in our case). We integrate the equations of motion defined by the Hamiltonians ([I^) 
and construct the surface, which consists of the gauge orbits, passing through to all the 
points of the evolution curves. The intersections of the surface with the level of the gauge 
fixing condition are just the evolution curves of the reduced system. Thus the explicit 
solutions of the evolution equations for the reduced system can be derived if the explicit 
description of the gauge orbits is known. In order to describe the gauge orbits we have to 
integrate the system of differential equations given by the kernel vector field (pUD 



-^Po = J23P2P3 -^x" = (x^p3 - X^P2] 



d d 1 

—Pi = = 

ds ds 

(17) 

—P2 = -poPs J-^X^ = -{POX^ + J23P3X°) 
-^P3 = P0P2 J^X^ = {PqX^ - J23P2X^), 

where s parameterizes the gauge orbit. Note that the r.h.s. of the differential equations 
for the momenta contains the momenta only. Hence we can integrate these equations 
separately. Then using the fact that the gauge transformations ( [Tl| ) are canonical ones 
we can resolve the rest equations for the coordinates with the help of the following formula 



dp^ 
dp. 



where Xq,p^ denote the initial conditions for the system (0). Note also that the phase 
coordinates x^, pi are not changed under the gauge transformations so in what follows 
we are going to omit them. It is easy to see that the system of differential equation for 
the momenta is similar to the system of equations for the ordinary SO3 top |1^. Using 



the analogy we can readily find two independent integrals of motion for the system (p!7| ) 

A = ^{P2Y + {P3? B = ^{p,Y + J{p,)\ (19) 

where J = J32 (in the rest of the paper we assume that J > 0). Let us consider the 
integrals ([T9|) as the new independent momenta. Then the rest third coordinate in the 
momentum space can be chosen in two different ways. First, one can introduce the 
following angle variable 

= arctg{ps/p2), (20) 
for which the momenta po;P2,P3 take the following form 



Pq = J B"^ — JA^sin?ip p2 = Acosif ps = Asimp (21) 
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Second, the third coordinate can be chosen as 

= arctg{y/jp3/po). (22) 

And in this case we get 



Po = Bcosip2 p2 = sjA^ - —sin?Lp2 P3 = -^sinip2 (23) 
The first way to choose the new momenta is apphcable when 

*^^<l. (24) 

because it is the case when the expression B^ — JA^sin^ip do not change the sign. Oth- 
erwise, one should use the second way of choosing coordinates in the momentum space. 
In both cases one can develop the projection method and find the evolution equations of 
the reduced system proceeding similarly at each stage. So in our paper we are going to 
discuss the projection method only for the first case ( |2^ ) omitting another possibility. In 
what follows, it is more convenient to use instead of Ob J Ob ^ 00 the coordinates, which are 
canonically conjugated to A, B and y?. They can be readily found in the following form 

2 n JAsin^if n 

xa = x cosip + x"" sirup ==x" 

V-D — J A'^sm^p) 

Bx^ 

'''' ^ ^B^ - JA^sm^^ ^^^^ 

2 • A JA^simpcosip r, 

x^ = -Ax smp + Ax^cosp t v> ■ ^ • 

In terms of A, B, p, xa,xb and x^ the solution for the system ([TtD looks as follows 
p{s) = am{Bs + F{po;k);k) A{s) = Aq B{s) = Bq 



_ / I - k'^sin^po 



xa{s) =x\ + ^^l- k'^sin'^ipo {Fk{p; k) - Fk{po; k))xl 

aVj, 



(26) 



Xb{s) =x%- sjl - k'^sin^ipQ {s + -^[Ffc(v9; k) - Fk{po] A;)])xJ 



where u = am{t] k) is the elliptic Jacobi function being the inverse function to the first 
order elliptic integral 

F{u- k) ^ 



'o Vl — k'^sin'^t 
Fk{u; k) denotes the corresponding derivative 

ksin'^tdt 



Fk{u; k) 



(1 — k'^sinHYl'^ 
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and AQ,Bo,ipo,x%x%x^ are the initial conditions of the respective phase coordinates. 
Now suppose we are given the initial conditions u^,vi,v?,V2 for the equations of motion 
of the reduced system ([l5| ) 



a 1 " a 1 

dt,u = — — dtX = 

Ot,U =V2 + 7 5 7-:t dt^U^ = JV2 

[ViW^ — JV2U^) 



^ti^^i = -7 ^ J TTi ^*2^i = 



(27) 



The initial point of the evolution curve (^) is mapped via the embedding ([13D to the point 
in the extended phase space with the phase coordinates x^,pu, such that x^ = pi = 
and the gauge fixing condition = is satisfied. Here we assume for simplicity that 
the inequality is also satisfied. The gauge fixing condition and the gauge invariant 
Hamiltonians (0) in terms of the phase coordinates (p!9|),(pOD,(25) take the following form 

AxACOsip = XipSimp, (28) 



C2 = Jl 



{Pl 



(29) 



2 2 ^ ^ 2 

Solving the equations of motion with respect to the Hamiltonians (^9]) we get that the 
momenta (p!9|,|20|) are conserved x^ and pi do not change their zero values and the coordi- 



nates (|25|) are linear functions of the times ti and t2, corresponding to the Hamiltonians 
Ci and C2 (|29|) respectively, namely 



Xn 











X 



^ x^ (po = ip 



xa + Atl Aq 



(30) 



A 



X 



B 



XB + Bt2 Bo = B. 



Here we suspend the phase space coordinates by and in order to show that they should 
be considered as the initial conditions for the system of differential equations (0) because 
in general a point of the phase curve (|30D does not satisfy the gauge fixing condition (|28 



but just define the orbit that should intersect the level of the gauge fixing condition in 
another point. Substituting now the phase space coordinates of the curve ( p6D with the 
initial data (BDf) into the gauge fixing condition (123) we get the following equation 



A{xA + Ah) 



X// 



1 — k'^sin'^if 
^ 1 — k'^sin'^(p{s) 
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-k^Jl - k^sin^ifi {Fk{(p{s); k) - Fk{(p; k)) 



(31) 



VIA 



where k = — = — . Since the r.h.s. of the equation (^) does not exphcitly depend on 
B 

the parameter s we can forget about it and treat (0) as the equation, defining (p as an 
imphcit function of ti. Really if we extract the dependence on cp in the r.h.s. of ( pT]) we 
get the function 



v^l — k'^sin'^ip 



kFk{ip; k) 



with the positive derivative 



> 0. 



Thus the intersection of the surface, which consists of the gauge orbits, passing through 
the points of the phase curve (|30|), with the level of the gauge fixing condition (p8|) is 
given by the following equations 



X^{tl, t2) = X 



1 — k'^sin'^if 



XA{ti,t2> = ? 



1 — kPsin'^ip 



XBih,t2) =Xb + Bt2 - ^^l-k^stn^^{F{ipihy, k) - F{^; k)) + 



(32) 



52 



-{Fk{ip{hy,k)-Fk{^;k)) 



A = A B = B (f = (f{ti), 



where (p{ti) is the implicit function given by (|3T|). Using the equations (0) we get the 
explicit solutions of the differential equations (pTj) in the following form 



52 - JA^sin'^ipih 



B 



u^{ti,t2) = XA{ti,t2) sirup (ti) + 



•^<fi (^1 ) ^ 



-XB(ti,t2 



JA 



^ COS ip{ti) + -g-XB(ti, t2)sinip{ti] 



(33) 



Vl(ti,t2 



B"^ — J A^sin^ipiti) ^2(^1,^2) = Asirvpiti), 
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where the initial data in the extended phase space can be obtained with the help of the 
following formulas 



A 



ip = —arctg- 



V 



- -2 - JAsin'^ip _i 

xa = u simp , _ u (34) 



B'^ — JA^siTn?p 



Xb 



B"^ — JA^sin^p 
J A^ sinpcos 

B"^ — jA'^sin'^p 



-r_2 - JA sinpcosp _i 
x<o = Au cosp u . 



4 Quantization 

In this section we consider the canonical quantization of the presented system restricting 
ourselves to the case of the positive J = J32 > , that is to the case when the quantum 
Hamiltonian H2 is an elliptic operator. We realize the Hamiltonians ([T5|) at the quantum 
level in such a way that the strong involution relation between them is not destroyed by 
quantum corrections. Then, given these commuting Hamiltonians, we solve the respective 
eigenvalue problem by requiring that the wave functions have no branch points. The 
surprising thing is that the spectrum of the Hamiltonian H2 turns out to be discrete and 
the coupling constant z/^ turns out to be negative. Our starting point is the classical 
Hamiltonians 



2J 2J{viu^ - V2v} 



Hi = -17T- + — TT^ ^2 = + (35) 



which are achieved from (|T^) by the obvious rescaling 



Vi Vi 



We perform the quantization of the system (|35|) in the momentum representation using 
the polar coordinates on the momentum plane 

v\ = rsincj) V2 = rcoscj). (36) 

In this setting 

■ , coscj) ^ sincj) 

ui = tn[stn(pOr H o^) U2 = in[cos(pOr o^) 
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■Oi = rsincp V2 = rcoscj), (37) 
and the quantum Hamiltonians take the following form 

2J 2JM2 '2' ^ ' 

where the operator M is defined as 

M = V1U2 — V2U1 = —ihdfj,. 

In order to determine the inverse operator to we use the following anzats for the basis 
wave functions 

^(r,0) = i?(r)<l>(0) (39) 
and represent the factor $(0) as the respective Fourier series 

00 



m=— 00 



2tt Jo 

Then the operator is naturally defined on the functions (^) by linearity if we set that 
the action of on the eigenfunctions e*™''^ of M looks as follows 

M-2^imc^ = ^—e*'"'^. (41) 

[nmy 



Of course we assume here that the zero mode in the expansion ( ^OD is vanishing $0 = . 
It is easy to see that thus defined Hamiltonians (^) are commuting operators and one 
may pose the respective eigenvalue problem 

H^^Ij = E^ij H21IJ = E2^. (42) 



Using the same anzats ( pQ]) for the wave function we see that this problem ( 142|) reduces 
to a single iteration equation for the coefficients 

m\^m+2 + - C<^m) + b^m = (43) 

since the dependence of the wave function (|3^) on r is obviously defined by the second 
equation in ( ^2]) 

R{r)=6{r-^2). 

The iteration equation (|43|) is equivalent to the following differential equation for the 
generating function $(z) = ^mz"^ 



1 



{zd^Yiiz' + -- cMz)) + b^z) = 0, (44) 

Z'^ 
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where 

In what follows it is more convenient to work with the unfolded form of this equation 
(^), namely 

z\z^ - cz^ + l)^"{z) + z{5z^ - cz^ - 3)$'(z) + {4z^ + 4 + = 0. (45) 



It goes with out saying that one can hardly find the solutions for the equation (^3|) in 
the explicit form. However, we may use the standard results of the complex analysis and 
conclude that the solutions are analytical functions on C besides may be the points where 
the polynomial P{z) = z'^{z'^ — cz^ + 1) is vanishing. Hence, we can qualitatively analyze 
asymptotical the behaviour of the solutions of ( ^51) in the vicinity of the roots of the 
polynomial -P(-2) = z^i^z'' — cz^ + 1) . And thus, we can derive the properties of the desired 
functions $(0), obtained by restricting the solutions to the unit circle {2 = 6*'^ : G R} . 
In fact we have to require that the restriction $((/>) is the function on the circle. It means 
that the respective analytic function $(2) have no branch points in the interior of the unit 
disc {2; G C : \z\ < 1} or the factors that arise from these branches should combine into 
the trivial factor 1 if the branch point is not the only one. In addition, we claim that the 
function $(0) belongs to L^(S'^) . The roots of the polynomial Pi^z^ = z'^{z^ — cz"^ + 1) that 
are worthy of notice are 2:1 = for arbitrary values of the parameters c and h 2:2,3 = ±1 
for c = 2 , 6 7^ 8 and 2:4,5 = ±2 for c = —2 and b ^ 8 . It is easy to show that the solutions 
in the vicinity of the point Zi = behaves like $(2) = z'^ , fi = 2 ± ^J^-b if 6 7^ —4 and 
behaves analytically if 6 = —4 . Hence 2:1 = is the branch point for our solutions if 
b 7^ —n^ , G N . It turns out that there are no other branch points in the interior of 
the unit disc and thus we conclude that the admissible values for b are —r? with natural 
numbers n . Considering the points ^2,3 = ±1 for c = 2 , 6 7^ 8 and 2:4,5 = ±z for c = —2 
and 6 7^ 8 it is easy to verify that the solutions of (|^) behaves like $(2;) = (2: — where 

w = 2:2,3 or w = 24,5 respectively and /? = . Since V n G N 



2 ^ 

we conclude that for an arbitrary b = —m? , n G N there exists a solution $(0) G L^{S^) . 
Note that in virtue of the relation 

b = TT = —n 

the coupling constant u"^ turns out to be negative since E2 , being the eigenvalue of the 
elliptic operator H2 is necessarily positive. Thus, we get the following spectrum for the 



Hamiltonians (|38|) 



2z/2 

El > E2 = n G N z/2 < 0. (46) 



The discrete spectrum of the Hamiltonian H2 is not so surprising because if we consider 
the exponential degeneration of our system ( |l6|) and quantize it in the representation of 
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the wave functions ip = 'ip{u'^, f i) we note that the momentum Vi defining the couphng 
constant of the two-particle Calogero-Moser system is quantized and thus the respective 
quantum Hamiltonian H2 also has a discrete spectrum. 



5 Concluding Remarks 

The main lesson we have learned from the above considerations is that we may not restrict 
ourselves to the Lie group symmetries in constructing the integrable systems via the 
symplectic reduction method. Of course if we go far beyond the Lie group symmetries the 
explicit description of the gauge orbits cannot be guaranteed and we cannot be sure that 
we are able to resolve the equations of motion via the projection method. Nevertheless we 
suspect that the non Lie group symmetries may be useful for constructing the mechanisms 
of the symplectic reduction for well known integrable systems. For example, although the 
mechanisms of the symplectic reduction for the Calogero-Moser systems associated to 



the root systems of the simple Lie algebras |l[] are known [|T^] the projection method for 
these system has not been developed. We suppose that the projection method for the 
systems can be presented if we construct the mechanisms of the symplectic reduction in 
the context of the non Lie group symmetries. It should be also very interesting to perform 
analogous symplectic reduction starting with the Odesskii bracket, which generalizes the 



Sklyanin bracket to the case of the arbitrary elliptic r-matrix . We suppose that 
the projection method for the respective systems should be related to solving the equations 
of motion for some higher dimensional tops ||13| as is the projection method for our system 
related to the equations of motion for the ordinary SO3 top. As it has been noted in the 
introduction the problem of finding the Lax representation for our system seems to be 
rather intriguing. Here we mention that this problem raises the question of whether it is 
possible to get the quadratic Sklyanin bracket via the reduction from the linear Poisson 
bracket. As far as we know the latter question is not answered yet even for the quadratic 



Poisson brackets, related to the constant triangular r- matrices |T6|. Finally, it is worth 
mentioning that in our paper we quantize the system starting with the reduced phase 
space of the model. However, it would be very interesting to obtain the same quantum 
mechanics via the respective quantum reduction, which may be performed with the help 
of the BRST approach [Q. We suspect that this quantum reduction should be some how 



related to the representations of the quantum Sklyanin algebra fl^] as are the quantum 
Calogero-Moser systems related to the representation of the respective Lie algebras 0]. 
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